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Abstract 

Two models of driven optical cavities, based on two-dimensional Ginzburg- 
Landau equations, are introduced. The models include loss, the Kerr nonlin- 
earity, diffraction in one transverse direction, and a combination of diffusion 
and dispersion in the other one (which is, actually, a temporal direction). Each 
model is driven either parametrically or directly by an external field. By means 
of direct simulations, stable completely localized pulses are found (in the di- 
rectly driven model, they are built on top of a nonzero flat background). These 
solitary pulses correspond to spatio-temporal solitons ( "light bullets" ) in the op- 
tical cavities. Basic results are presented in a compact form as stability regions 
in a full three-dimensional parameter space of either model. The stability region 
is bounded by two surfaces; beyond the left one, any two-dimensional (2D) pulse 
decays to zero, while quasi-lD pulses, representing spatial solitons in the optical 
cavity, are found beyond the right boundary. The spatial solitons are found to 
be stable both inside the stability region of the 2D pulses (hence, bistability 
takes place in this region in the two models) and beyond the right boundary 
of this region (although they are not stable everywhere). Unlike the spatial 
solitons, their quasi-lD counterparts in the form of purely temporal solitons are 
always subject to modulational instability, which splits them into an array of 
2D pulses, that further coalesce into two final pulses. A uniform nonzero state 
in the parametrically driven model is also modulationally unstable, which leads 
to formation of many 2D pulses that subsequently merge into few ones. 



1 Introduction 

Complex Ginzburg-Landau (CGL) equations is a class of universal models de- 
scribing pattern formation in media of various physical nature, combining non- 
linearity with dissipative and dispersive linear properties. They find important 
applications to traveling- wave convection (see Ref. |lj and references therein), 
fiber optics M, and in other fields. Nonlinear terms in a CGL equation, as 
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well the linear ones, may be both dissipative, accounting for loss and gain, and 
conservative. 

An important class of patterns that appear in different versions of the CGL 
equations are solitary pulses (SPs). In particular, a GL model with the cubic- 
quintic (CQ) nonlinearity was first introduced by Petviashvili and Sergeev Q 
in the two-dimensional (2D) case, with an objective to construct stable fully 
localized 2D pulses. The nonlinearity of the CQ type was necessary, as the 
stability of a pulse requires, first of all, the zero background to be stable, which 
implies that the linear part of the equation could not contain gain, hence a gain 
term had to be placed in the cubic nonlinearity. Finally, to provide for the 
overall stability of the model, the cubic gain had to be capped with a quintic 
lossy term. 

Stable SPs ("solitons") in 2D CGL equations of various types were found 
in Refs. Q. Stable axisymmetric solitons with an internal vorticity ("spin") 
have been found too in the isotropic CQ CGL equation || . On the other hand, 
the 2D CGL equation may find its most plausible experimental realization as a 
model governing evolution of optical spatiotemporal pulses in a planar nonlinear 
waveguide (2D optical cavity); in that case, however, the equation is strongly 
anisotropic due to the natural difference between the propagation direction (z) 
and the transverse one (x) along which the electromagnetic wave diffracts. In a 
rather general form, the corresponding anisotropic 2D model was put forward 
inRef. §: 

1 i 

iu z + -u xx - (- + P)u tt + [iu + (1 - «7i) \u\ 2 u + ij 2 \u\ 4 u] = 0. (1) 

In this equation, u(z,x,t) is the local amplitude of the electromagnetic field, 
t = T — z/Vq is the so-called reduced time, where T is time proper, and Vo is 
the group velocity of the carrier wave propagating in the z-direction fej. Further, 
coefficients in front of the terms accounting for the transverse diffraction (^u xx ), 
linear loss (iu), optical filtering (— hiiitt, which formally looks like diffusion 
in the t-direction), and Kerr nonlinearity (|m| 2 m) are all normalized to be 1. 
The remaining free parameters (3, 71, and 72 control the temporal dispersion, 
nonlinear cubic gain, and quintic loss (71 and 72 must be positive, while (3 may 
have either sign or be equal to zero). SP solutions to Eq. (|l|) in the form 
u(z, x, t) = exp(jfcz) v(x, t), with the function v(x, t) exponentially vanishing at 
\x\ — > 00 and \t\ — > 00, were found in Ref. ||, and their stability region was 
identified in the corresponding parameter space. 

The model based on Eq. ([j]) assumes that the planar waveguide carrying 
the field is equipped with uniformly distributed linear bandwidth-limited gain 
(which induces the filtering term) combined with a saturable absorber that 
together give rise to the CQ nonlinearity. For applications to nonlinear optics, it 
is also quite interesting to consider models where wave patterns are supported, 
instead of the intrinsic gain, by an external pump field, that leads to driven 
CGL equations, in which it is sufficient to take into regard the ordinary cubic 
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(Kerr) nonlinearity, while the linear diffraction, dispersion, and filtering terms 
remain the same as in Eq. (|l]). In fact, two different driving terms may appear 
in the accordingly modified equation: a parametric drive and direct forcing, 
see below. It is quite interesting to find 2D spatiotemporal SPs in the driven 
optical models, as they have a real chance to be observed in the experiment, 
and are also of interest in their own right as completely localized pulses in the 
anisotropic 2D CGL equations of the new type. 

In this work, we find these pulses and demonstrate their stability. Besides 
that, we also study quasi-lD solitons of two different types, which represent spa- 
tial and temporal solitons in the underlying optical models. It will be demon- 
strated that the spatial solitons are stable in a broad parametric region (although 
not everywhere); in particular, they coexist as stable attractors with the stable 
2D pulses, hence both models are bistable in the corresponding regions. On 
the contrary to that, the temporal solitons are always modulationally unstable, 
and simulations show that the instability splits them into an array of 2D pulses, 
which subsequently coalesce into two 2D solitons. A uniform nonzero state 
in the parametrically driven model is also subject to modulational instability, 
which leads to a formation of a large number of 2D SPs; later, they merge into 
few final 2D solitons. 

The rest of the paper is organized as follows. The two versions of the driven 
model are introduced in section 2, and SP solutions found in them are displayed 
in section 3. The results of extensive numerical simulations are summarized in 
the form of explicit stability diagrams in the corresponding three-dimensional 
parameter spaces. The paper is concluded by section 4. 



2 The driven models 

2.1 The parametric drive 

In the ID case, the parametrically driven CGL equation takes the form which 
generalizes the earlier studied parametrically driven damped nonlinear Schrodinger 
(NLS) equation (see Ref. 0, in which exact SP solutions to the latter equa- 
tion were found, and their stability was investigated, and further results in 
Refs. H|). The generalization to the 2D case of the same type as in Eq. (|l]), 
which describes a planar waveguide (optical cavity) with the Kerr nonlinearity, 
is straightforward: 

\ % 1 

iu z + -u xx - (- + (3)u u + -j\u\ 2 u = (k - i)u + ju*, (2) 

where the asterisk stands for the complex conjugation. The gain, which is 
accounted for by the last term in Eq. (j^) with a real coefficient 7 is provided by 
a pump wave copropagating with the u-field at the double frequency. The real 
coefficient k determines a wavenumber mismatch between the signal and pump 
waves. 
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Dropping the filtering term and assuming < (which physically corre- 
sponds to anomalous dispersion of the electromagnetic waves in the waveguide 
0) transform Eq. (||) into a 2D isotropic damped parametrically driven NLS 
equation, which is used as a phenomenological model of the Faraday ripples in 
hydrodynamics |9| , in which case the evolutional variable z in Eq. (|^) is re- 
placed by time T, and the former temporal variable t is replaced by the second 
transverse coordinate y. A more general isotropic equation, which includes the 
diffusion term, 

iu z + (~ - iD)\7 2 u + ^\u\ 2 u = (k - i)u + ju* , D > 0, (3) 

is also used as a model to describe localized 2D states, the so-called oscillons, 
observed in experiments with the Faraday resonance 0. Getting back to the 
damped isotropic NLS equation with the parametric drive, we note that this 
equation (with the evolutional variable z) finds applications to the description 



of spatial patterns in 3D driven optical cavities 1 1 . In a recent preprint p2| ] , it 
has been shown that axisymmetric solitons may be stable in the latter model, 
despite the presence of the collapse in its undamped undriven version (which is 
the NLS equation proper). 

In the present work, the objective is to find stable SPs in the strongly 
anisotropic full CGL model (||), which will have the meaning of spatiotemporal 
optical solitons [|l3) propagating in the planar waveguide. It seems obvious that 
pulses may only exist if the gain parameter exceeds the loss coefficient, i.e., 
7 > 1. On the other hand, an obvious necessary condition for the stability of 
pulses is the stability of the zero background against perturbations of the form 

u = aexp (az + iqx — iu)t) + 6exp [a* z — iqx + iuit) (4) 

with infinitesimal amplitudes a and 6, where q and to are arbitrary real wavenum- 
ber and frequency of the perturbation. A dispersion relation which determines 
the (complex) instability growth rate a as a function of q and uj can be easily 
found from the linearized equation (0): 



l -q 2 ~pu 2 + k 
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Re a + 1 + -uj 2 I + ilma 



r, (5) 



which gives rise to two different branches of the dispersion relation, viz., Re a = 
— (l + ^w 2 ) with Imcr ^ 0, and the one with Imcr = and a determined by 
the equation 

+ 1 + ^ 2 ) 2 ' = l 2 -(^q 2 -^ 2 + kj \ (6) 

Obviously, the former branch always satisfies the stability condition, Re a(q, ui) < 
0. Straightforward analysis of Eq. (p|), with regard to the above assumption 
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7 > 1, shows that, if uj = 0, the stability condition for the latter branch amounts 
to the well-known inequalities M, || 7 2 < 1 + fc 2 , fc > 0. 

The consideration of Eq. rtq) with u demonstrates that, in this case, 
the most dangerous perturbations are those with q — 0. If (3 < 0, they do not 
generate any additional instability, and if (3 is positive, a new stability condition 
appears, 7 2 < (fc + 2/3) 2 /(l + 4/3 2 ); it is easy to check that the above-mentioned 
necessary condition 7 2 < 1 + fc 2 is a straightforward corollary of this inequality. 
Thus, the eventual system of the conditions providing for the stability of the 
zero solution to Eq. (||) (provided that the system is infinitely extended in both 
x— and i-directions) is fc > and 

l 2 <l + k 2 , if/3<0; 
7 2 < (fc + 2/?) 2 /(l + 4/3 2 ), if/3>0. [l> 

It is also relevant to consider nonzero spatially uniform solutions to Eq. ( 
||), which have the form 

u = A Q exp (i(f>) , (f> = --sm- 1 (l/j), A = ^2 [fc + jcos(2(j>)}. (8) 

This uniform solution exists for 7 > 1. The consideration of the modulational 
stability of the solution against small perturbations with the wavenumber q and 
instability growth rate a (cf. Eq. (|])) yields a dispersion relation 

2 1 1 \ 

= 7 2 -7^oCOs(2</>)+^/4. 



Qq 2 - (3uj 2 + k - Alj + (rcct + 1 + io> 2 ^ + ilmcr 



As well as Eq. (||), this equation has a solution branch with Imcr = 0, then 

a = -1 - uo 2 /2 ± - 7A 2 cos(2</>) + A%l 4 - (q 2 /2 - ptu 2 + k - A 2 ) 2 . (9) 

The maximum value of this expression, which is attained at u> = and q = 
\/2(k — Aq), is always positive, cr max = — 1 + \/l + fc 2 - Therefore, the uniform 
state (^|) is always unstable against modulations in the x-direction, the same 
way as in the ID damped parametrically driven NLS equation |). 



2.2 The directly forced model 

The light signal in a lossy waveguide can also be directly (rather than paramet- 
rically) supported by a pump wave launched at the same (rather than double) 
frequency. For spatial patterns in 3D lossy optical cavities (a typical example 
is a domain wall |l4| ) , this gives rise to the well-known isotropic directly forced 
2D CGL equation, with the propagation coordinate z playing the role of the 
evolutional variable, see a review pi (a similar model describes a cavity with 
a trapped Bose-Einstein condensate in the presence of source and sink of atoms 
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[[16)). The optical patterns generated this way may find an application in the 
design of optical memory 0, although fully localized pulses may be subject to 
collapse under the action of the self-focusing Kerr nonlinearity pl| . 

The governing equation for the spatiotemporal optical signal in the directly 
forced nonlinear waveguide differs from Eq. (Q) by the driving term: 

1 i 1 

iug + -u xx - (- + f3)u tt + -\u\ 2 u = (k- i)u - /, (10) 

where / is the amplitude of the forcing field. Obviously, this equation can 
only produce localized patterns on top of a constant background u = uq, 
whose amplitude is determined by a cubic equation following from Eq. (|l0|), 
(1/2) \uq\ 2 uq + (i — k) uq + f = 0. The background can be subject to modula- 
tional instability ^ ; however, if the amplitude uo is small enough, the instability 
will be so weak that it may be disregarded for finite values of the propagation 
distance that are relevant to the experiment. 



3 Numerical results 

3.1 Two-dimensional pulses in the parametrically driven 
model 



Equations (|2j) and (10) were solved numerically by means of a pseudospectral 
code with 256 x 256 modes and periodic boundary conditions in x and t, with 
a fixed size of the system in both directions, L x = L t — 20. As well as in Ref. 
[||, a stable 2D SP was generated for the first time, starting with an initial 
Gaussian pulse localized in both x and t, with its center placed at the center 
of the integration domain. In order to generate stability diagrams (see below), 
we then varied the parameters by small steps, the initial configuration for each 
simulation being the stable pulse produced by the simulation at the previous 
step. 

In Fig. 1 we display a typical example of a stable 2D completely localized 
pulse, which was produced by the numerical solution of Eq. (||), following the 
procedure outlined above. Results of extensive simulations, which were carried 
out for many different values of the parameters, can be presented in a compact 
form as a stability diagram for the 2D pulses in the three-dimensional parametric 
space (7, fc, /?), which is displayed in Fig. 2. We stress that this parameter space 
comprises all the coefficients of Eq. (^). 

The stability region proper is located under the "roof" in Fig. 2. To the 
left of the stability border which is shown by crosses, no localized 2D pulses can 
be created: in this case, an initial 2D pulse configuration quickly decays into 
the stable zero state. On the other side, i.e., to the right of the border denoted 
by rhombuses, localized 2D pulses were not found either. However, the result 
of the evolution is different in this case: the initial 2D pulse expands along the 
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i-axis to generate a quasi-lD pulse, which is localized in the x-direction, and 
delocalizcd in the t -direction. Such a quasi- ID solution to Eq. (||) can be found 
in an exact form, 

u(x,t) = VZne^sechinx), (ft = ~ sin _1 (l/7), k= \/2 [k + 7Cos(2c/>)] . (11) 



In terms of the underlying optical models, the solution (11) corresponds 
to a spatial soliton. A typical example of such a soliton is shown in Fig. 3. 
Straightforward comparison shows that the numerically found quasi- ID pulse 
completely coincides with the analytical solution (|ll]) . Simulations demonstrate 
that this spatial soliton is stable everywhere in the stability region of the 2D 
pulses shown in Fig. 2, i.e., the model is bistable in this region. Beyond the 
right stability border of the 2D pulses, shown by rhombuses in Fig. 2, i.e., in 
the region where only the spatial solitons are found, they remain stable. 

For comparison with these results, we have also performed direct simulations 
of the isotropic equation (||). Stable 2D localized pulses exist in this model too, 
Fig. 4 displaying an example for D — 0.1. For D — 0.1 and k — 1.3, the stable 
axially symmetric 2D pulses have been found in the interval 1.20 < 7 < 1.57. 
For 7 < 1.20, the initial 2D pulse decays to the zero state, while for 7 > 1.57 the 
initial 2D pulse evolves to the spatially uniform state (||). In this connection, 
we note that a critical branch of the modulational-instability growth rate for 
this uniform solution to Eq. @ can be found, as a function of the perturbation 
wavenumber q, in the form a(q) = — 1 — Dq 2 + y/l + k 2 — (q 2 /2 + k — A0) 2 , cf. 
Eq. (||). It is easy to check that, for D = 0.1 and k = 1.3, a(q) is negative for 
any q, hence the uniform state is stable. 

3.2 Modulational instability of quasi-one-dimensional pulses 
in the parametrically driven model 

In the context of the results presented above, the spatial solitons (quasi-lD 
pulses) were stable. However, they become unstable against t -dependent per- 
turbations for sufficiently large negative values of the dispersion coefficient (3 in 
Eq. (||) , which corresponds to anomalous dispersion in terms of nonlinear optics 
(in fact, modulational instability of spatial solitons against temporal perturba- 
tions in a system with anomalous dispersion is a general effect which has been 
known for a long time fl9|| ). In particular, for k = 1.3 and 7 = 1.3, this insta- 
bility occurs if j3 < —1.08. Figure 5 displays a result of numerical simulations 
in this case. The initial condition including a small perturbation is taken as 



u(x, t,z = 0) = V2h 



/8ttA /27rfV. 
1 + 0.05 cos I — 1 + 0.01 cos I — 1 ) cos (ft + i sin < 



sech 



(12) 

where K and (ft are the same constants as in Eq. ([ll]). Figure 5(a) displays the 
evolution of \u(x = L/2, t)\ at j3 = —2.5. While the pattern remains localized in 
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the x-direction during the evolution, it is obvious in Fig. 5(b) that the growth 
of the modulation splits the original spatial soliton into four 2D spatiotemporal 
pulses, but three of them merge together, thus two 2D pulses survive (splitting 
of a spatial solitary beam into temporally localized pulses under the action of 
modulational instability was earlier observed in numerical (2^, and labora- 
tory |2l[| experiments with second-harmonic-generating optical systems) . Figure 
5(b) displays an eventual established state including the two 2D pulses. This 
modulational instability is a generic type of transverse instability for quasi-lD 
pulses (another generic type is the zigzag instability found in Ref. || in the 
context of the 2D GL equation with the CQ nonlinearity) . 

Another type of quasi- ID pulse solution to Eq. (||), which is localized in 
the ^-direction and uniform in the x-direction, is also possible. Such an x— 
independent solution, u(z,t,z) = U{z,t), represents a temporal soliton in the 
optical cavity, and satisfies the equation 

%U Z -C-+ f3)Utt + \\U\ 2 U = {k- i)U + ~fU* (13) 

(an exact analytical solution to this equation is not available, unless f3 = 0). 
For the same case, k — 1.3 and 7 = 1.3, as considered above, this type of the 
quasi-lD pulse solution exists and is stable against i-dependent perturbations if 
(3 < —0.812. Figure 6(a) displays such a pulse at (3 = —1, k = 1.3 and 7 = 1.3, 
with its center at the point t = L/2. To study its stability against ^-dependent 
perturbations, we took an initial condition 

u(x, t,z = 0) = [l + 0.05 cos(27ra/L)] Re U{t) + ilm U(t), 

where U(t) is the unperturbed z -independent SP solution to Eq. ([l3] ). Mod- 
ulation in the x- direction grows and splits the temporal soliton into many 2D 
spatiotemporal pulses, as it is evident in Fig. 6(b), which displays a 3D plot of 
|it(x, £, z)| at z = 200. Such a temporal soliton, extended along the x-direction, 
seems to be always unstable against the x-dependent perturbations, breaking 
up into a large number of 2D pulses. 

One can also look for more general quasi-lD pulses of the form u(x,t,z) = 
U(x — ct), which may be interpreted as either moving spatial solitons, or oblique 
temporal ones, the function U satisfying an equation 

(1 - 2c 2 p - ic 2 ) U" + \U\ 2 U = 2(k - i)U + 2 7 C7*. 

A detailed study of such general quasi- ID solitons is beyond the scope of this 
work, but one may assume that there is a critical value of the velocity c sepa- 
rating completely unstable solutions and those which may be stable. 

We have also simulated the development of the modulational instability of 
the spatially uniform state (H) in the parametrically driven model (||) (the pres- 
ence of this instability was demonstrated analytically in the previous section). 
Formation of many 2D pulses may be expected, at sufficiently small /3, as a 
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result of the instability development. As an initial condition, we took the uni- 
form state (^) randomly perturbed by an initial small disturbance. Figure 7 
displays a snapshot of the pattern obtained at z = 600 for (3 — —2.5, k = 1.3 
and 7 = 1.3. At first, many 2D pulses are created; however, they subsequently 
merge, so that only a few 2D pulses eventually survive, which is obvious in Fig. 
7. 

3.3 The directly forced model 

Numerical simulations of the directly driven CGL equation ([h]) also readily 
produce stable 2D pulses (built on top of the small-amplitude background, see 
above). A typical example of such a pulse is shown in Fig. 8. As well as 
in the parametrically driven model, in this case extensive simulations of Eq. 
( fl0| ) have made it possible to identify a stability region for 2D SPs in the full 
three-dimensional parameter space of the model, (/, k, (3). This region is located 
under the "roof" in the parameter space, see Fig. 9. Note that the parameter 
space in Fig. 9, as well as that in Fig. 2, incorporates all the parameters of the 
model. 

The two boundaries of the stability region in Fig. 9 have the meaning similar 
to that of the stability boundaries shown in Fig. 2 for the parametrically driven 
model: to the left of the boundary shown by crosses, any initial pulse decays 
to zero, and to the right of the second boundary, shown by rhombuses, the 
fully localized 2D pulse evolves to a quasi- ID pulse oriented along the f-axis 
(the spatial soliton, in terms of the underlying optical models). As well as in 
the parametrically driven model, in the directly forced one the spatial solitons 
coexist as stable patterns with the 2D pulses (hence, the model is bistable in 
that region), and they remain stable across the right border of the stability 
region of the 2D solitons, i.e., across the surface shown by rhombuses in Fig. 9. 



4 Conclusion 

In this work, we have introduced two different two-dimensional models of the 
Ginzburg-Landau type with the Kerr (cubic) nonlincarity. The corresponding 
equations are driven either parametrically or directly by an external field. The 
models describe spatio-temporal dynamics in nonlinear dispersive lossy optical 
cavities, where the loss also includes a filtering term. Both models are strongly 
anisotropic, featuring only diffraction in one (spatial) direction, and a combi- 
nation of effective diffusion and dispersion in the other (temporal) direction. 
By means of direct simulations, we have found stable two-dimensional solitary 
pulses in both models (in the directly driven one, the pulses are built on top of 
a nonzero flat background). These two-dimensional solitary pulses correspond 
to spatio-temporal solitons ("light bullets") supported by the driving field in 
the planar nonlinear waveguide. 
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The main results were presented in a compact form as stability regions in the 
full three-dimensional parameter space of either model. The stability region is 
bounded by two surfaces. Beyond one of them, any initial pulse decays to zero, 
while stable quasi-one-dimensional pulses, extending along the temporal direc- 
tion (i.e., spatial solitons, in terms of the optical cavity), were found beyond the 
other boundary. The spatial solitons are also stable inside the stability regions 
of the two-dimensional pulses, so that both models demonstrate bistability in 
these regions. Unlike the spatial solitons, their quasi-one-dimensional counter- 
parts in the form of temporal solitons are always modulationally unstable. It 
was demonstrated that the instability splits temporal solitons into an array of 
spatiotemporal two-dimensional pulses, which, in turn, coalesce into two final 
pulses. Similarly, the instability of the uniform nonzero state in the paramet- 
rically driven model leads to formation of many two-dimensional pulses, which 
subsequently merge into few two-dimensional solitons. 
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Figure Captions 

Fig. 1. An example of the stable stationary two-dimensional solitary pulse, 
found as a solution of the parametrically driven complex Ginzburg-Landau 
equation (Q). The field |u| is shown vs. the temporal and transverse spatial 
coordinates t and x. Values of the parameters are k = 1.3, (3 = —0.3, 7 = 1.30. 

Fig. 2. The region of stable two-dimensional solitary pulses in the full pa- 
rameter space of Eq. (@) . The stability region is located under the "roof" . To 
the left of stability region, only zero solution is possible, and to the right of it, 
the 2D pulse evolves to a stable quasi-one-dimensional pattern, see an example 
in Fig. 3. 

Fig. 3. An example of a stable quasi-one-dimensional pulse found in the para- 
metrically driven complex Ginzburg-Landau equation (^) . Values of the param- 
eters are k = 1.3, j3 = —0.3,7 = 1-36 (a point corresponding to these values is 
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located slightly to the right of the stability region shown in Fig. 2). 

Fig. 4 An example of a stable axisymmetric two-dimensional solitary pulse 
found as a solution of the parametrically driven isotropic complex Ginzburg- 
Landau equation (||). The field \u\ is shown vs. the temporal and transverse 
spatial coordinates. Values of the parameters are k = 1.3, D = 0.1 and 7 = 1.30. 

Fig. 5 (a) Evolution of \u(L/2,t, z)\ for the parametrically driven complex 
Ginzburg-Landau equation (g). The initial state is a ID localized pattern ( [l2| ) 
(spatial soliton) including a small modulational perturbation, values of the pa- 
rameters being k = 1.3,(3 = —2.5,7 = 1-3 ■ (b) The established stationary 
pattern including two 2D pulses. 

Fig. 6 (a) The quasi-one-dimensional pulse, which is localized in the t -direction 
around t = L/2 (temporal soliton), at /3 — — 1, k = 1.3 and 7 = 1.3, in the para- 
metrically driven complex Ginzburg-Landau equation ( ^) . (b) The result of the 
instability development of this soliton, shown as a plot of |u(a;, t)| at z = 200. 

Fig. 7 Patterns produced by the modulational instability of the uniform state 
in the parametrically driven complex Ginzburg-Landau equation (Q) at z = 600. 
Values of parameter are (3 = —2.5, k = 1.3 and 7 = 1.3. 

Fig. 8. An example of the stable stationary two-dimensional solitary pulse, 
found on top of the nonzero flat background, in the directly forced complex 
Ginzburg-Landau equation (|l0|). The field \u\ is shown vs. the temporal and 
transverse spatial coordinates t and x. Values of the parameters are k = 3, (3 = 
-0.3,/ = 3.1. 

Fig. 9. The region of stable two-dimensional solitary pulses (built on top of 
the nonzero flat background) in the full parameter space of the directly forced 
complex Ginzburg-Landau equation (|lo|). The stability region is located un- 
der the "roof". To the left of stability region, the 2D pulse decays to the 
flat-background solution, and to the right of the stability region, the 2D pulse 
evolves to a stable quasi-one-dimensional pattern (spatial soliton). 
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